Shell model test of quadrupole properties predicted by the 
rotational formula — the degenerate SDI interaction and 

non-degenerate FPD6 

A. Escuderos and L. Zamick 
Department of Physics and Astronomy, 
Rutgers University, Piscataway, New Jersey 08854 
(Dated: February 9, 2008) 

Abstract 

In the rotational model for a K = band in an even-even nucleus, there is a single parameter — 
Qo, the intrinsic quadrupole moment. All S(£J2)'s in the band and all static quadrupole moments 
are expressed in terms of this one parameter. In shell- model calculations, this does not have to 
be the case. In this work, we consider ground-state bands in 44 Ti, 46 Ti, 48 Ti, 48 Cr, and 50 Cr 
with two different models. First, we use a Surface Delta Interaction with degenerate single-particle 
energies (SDI-deg). We compare this with results of a shell-model calculation using the standard 
interaction FPD6 and include the single-particle energy splitting. Neither model yields a perfect 
rotational 1(1 + 1) spectrum, although the SDI-deg model comes somewhat closer. Overall, the 
simple rotational formula for B(E2)'s and static quadrupole moments hangs together very nicely. 

PACS numbers: 21.60.Cs,21.60.Ev,27.40.+z 
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I. INTRODUCTION 

In a previous publication [lj], we considered the relationship between the static quadrupole 
moments of the 2± states and the corresponding B(E2)'s in the sd and fp shells. We used the 
shell model to operationally define a ratio of intrinsic quadrupole moments Q (S)/Q (B), 
where Qo{S) is obtained from the static quadrupole moments and Qq(B), from the B(E2) , s. 
In the simple rotational model, this ratio would equal unity, and in the harmonic vibrator 
model, it would equal zero. 



There have been tests o 



model by Bender et al. 0, ^j]. 
quadrupole deformations by Jaqaman jmd Zamick 
and more recently by Sagawa et al 



this ratio in other models such as the Skyrme Hartree-Fock 
This same model has been used to test the systematics of 



4|, Zheng et al. [5|, Retamosa et al. p, |7[, 
At the same time as Ref. Ill], there appeared an 
article on the same topic but with a different approach by S.M. Lenzi et al. |9J and, more 



101 ] . Very recent references 



111 . [12I show ever increasing 



recently, by G. Thiamova et al. 
interest in this subject. 

In another vein, random interaction studies were performed by Velazquez et al. [13( and 
by Zelevinsky and Volya [14J]. They found two spikes (i.e., high probabilities) in the Alaga 
ratio A = 5Q 2 /[16irB(E2)o->2] at A = and A = A = 4/49, which can be associated with 
the vibrational and rotational limits, respectively. 

On the other hand, for the nuclei that we considered in Ref. [lj], the experimental ratio 
Qo(S) / Qo(B) was, for the most part, large, sometimes exceeding one, e.g., for 20 Ne and 
50 Cr. The one exception was 40 Ar, where the ratio was 0.06. 



151 ] have shown that, in a full fp space, one 



We would like to mention that Poves et al. 
gets a quasi-rotational band in 48 Cr. The spectrum is not /(/ + 1) exactly but there are 
strong intraband E2 transitions between the levels. 

In this work, we shall extend this study by considering states of higher angular momentum 
as well. All calculations have been done using the shell-model code ANTOINE jlfij ]. 

We define 

RsB = q°(g>" = 1 q^) = — Q< 2+ > , (1) 

<5o(B)2-o 2yiftF v /B(£'2) M) V s ( B2 )2-o 

and 

M(O) _Qo(S)i _22I + 3Q(I) 



where Q(I) is the (laboratory) static quadrupole moment of a state of angular momentum 

I. Note that the Alaga ratio mentioned above is A = -^\Rsb(^)\ 2 ■ 
We also define 

M(B),^ 2 = Q ^ ] '-- 2 • (3) 
Now for a K = rotational band, we have 

B(E2) I ^ I _ 2 = -|-(J200|J - 2, 0) 2 Q|(/) , (4) 

107T 

Hence, 

Mm] 2 2 (2/ -1)(2/ + 1)^(^2)^-2 

M(5)m " 15 1(1-1) B(E2) 2 ^ ■ (6) 

Note that both M(Q) 2 and M(B) 2 ^o are equal to 1 by definition. 

II. PRELIMINARY REMARKS 

A. The Schematic Surface Delta Interaction 

In the past, schematic models, despite giving somewhat oversimplified descriptions of the 
structures of nuclei, proved to be invaluable in casting insights into the trends of nuclear 
structure. As an example, Elliott's SU(3) model showed how one could approach rotational- 



model spectra in the Shell Model [17j. This was shown with a two-body momentum- 
dependant long-range quadrupole-quadrupole interaction and did not include the effects of 
spin-orbit splitting. Indeed, SU(3) models emanating from this interaction are emphasized 
in Refs. QQ. 

In this work, the main thrust will be to use a realistic interaction with correct single 
particle splittings. However, since our results for R$b agree with the rotational model even 
when the spectra are not rotational, we are motivated to get insight into this result by using 
a different schematic interaction, one which does not have 1(1 + 1) spectra but still exhibit 

n 

collective behaviour. In contrast to Elliott's long range interaction [171 ] . we will use the 
surface delta interaction of Moszkowski with degenerate single particle energies [l^ |. 

In Fig. [1], we show the calculated spectrum of 48 Ti with a surface delta interaction as 
compared with the /(/+!) rotational spectrum. The parameters have been adjusted so that 
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the excitation energies of the first 2 + states are the same. The rotational spectrum is more 
spread out but the surface delta spectrum still has rotational features and the spectrum is 
actually closer to the truth than the /(/ + 1), at least for 48 Ti. Hence, the surface delta 
interaction will serve as a good counterpoint to the more realistic interaction considered in 
the next section. 
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FIG. 1: Comparison of the SDI-deg and rotational 7(7 + 1) spectra for Ti. Both interactions 
have been fitted to reproduce the experimental splitting of the Of and 2f states. 
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B. Experimental Comparison 

In Ref. S.M. Lenzi et al. compare the B(E2)'b of N = Z nuclei, 44 Ti and 48 Cr, 
with the SU(3) model of Elliott making use of the momentum-dependent quadrupole- 
quadrupole interaction. Although the SU(3) model gives a rotational spectrum, the S£ ! (2)'s 
do not follow the rotational formula. Rather, with increasing J, they decrease relative to 
the rotational formula. The authors make the point that, despite the absence of spin-orbit 
splitting, the SU(3) results are not too bad compared with experiment. For 48 Cr they refer 
to the experiments of Brandolini et al. [19]. We shall also use this reference for our analysis 
of Qo(B) I ^ I ^/Q (B) 2 ^ . 



In tables HI and ITU we can see the results from Ref. 19j for 48 Cr and 50 Cr, respectively. 
Note that, in the case of 48 Cr, M(i?)/_ > /_2 decreases almost linearly with /. There is more 
experimental information, but, beyond the results shown in the tables, things get lost in 
band crossing. 



TABLE I: Experimental B(E2)'s [e fm ] and ratio M(B)i^j-2 (see text) for Cr; experimental 



data are taken from Ref. 



191 ]. except the B(E2)2-+o, which is taken from Ref. [2C] 







B(E2) exp [e 2 fm 4 ] 


M(fl)/_>/_ 2 


2 - 


♦ 


272 


1.00000 


4 - 


♦ 2 


329(110) 


0.92016 


6 - 


♦ 4 


301(78) 


0.83863 


8 - 


♦ 6 


230(69) 


0.71651 


10 - 


8 


195(54) 


0.65098 


12 - 


♦ 10 


167(25) 


0.59716 


14- 


♦ 12 


105(18) 


0.47057 


16 - 


♦ 14 


37(8) 


0.27805 



III. RESULTS 



As in Ref. [l|, in Table [TTT1 we check the relationship between the B(E2) 2 +__ rQ + and the 
static quadrupole moment (3(2+). 



TABLE II: Same as Table U but for 50 Cr. 

B(E2) cxp [e 2 fm 4 ] M(B)j^ 2 

2^ 216 1.00000 

4^ 2 204(57) 0.81309 

6^4 235(47) 0.83154 

8^6 205(51) 0.75909 

10(1) -> 8 72(14) 0.44389 

10(2) -> 8 131(26) 0.59875 

TABLE III: Q(2+) [e fm 2 ], B(E2) 2 +^ + [e 2 fm 4 ], and R SB for 44 > 46 > 48 Ti and 48 < 50 Cr obtained from 
a full fp shell calculation with the interactions SDI and FPD6. 







SDI-deg 






FPD6 






g(2 + ) 


B(E2) I 


Rsb 


Q(2+) 


B(E2) I 


Rsb 


44 Ti 


-26.319 


165.74 


1.0092 


-20.156 


121.45 


0.9029 


46 Xi 


-29.349 


206.06 


1.0093 


-22.071 


136.41 


0.9329 


48 Xi 


-32.337 


247.38 


1.0149 


-17.714 


112.16 


0.8257 


48 Cr 


-37.936 


377.60 


0.9638 


-33.271 


275.68 


0.9892 


50 Cr 


-41.681 


435.82 


0.9856 


-30.955 


243.80 


0.9787 



In the rotational model, 



Q(2+) = -lQo(S) , (7) 
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Qo(Bf 



B ^t-o t = ^ ■ (8) 
Furthermore, the ratio Rsb = Qo(S)/Qo(B) should equal 1. 

In Table HTT1 we give Rsb for two models: a surface delta interaction with degenerate single- 
particle energies (SDI-deg) and the more realistic calculation with the FPD6 interaction 
including single-particle energy splittings, both in a full fp space. The strength of the SDI- 
deg interaction was chosen to fit the experimental excitation energy of the 2f state. The 
nuclei considered are 44 Ti, 46 Ti, 48 Ti, 48 Cr, and 50 Cr. Note that, from Figs.[2H6] neither SDI- 



6 



deg or FPD6 have pure rotational spectra, although SDI is closer to a rotational spectrum, 
undoubtedly due to the fact that there are no single-particle splittings. 
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FIG. 2: Full /p space calculations of even-J states in 44 Ti with the SDI-deg and FPD6 interactions, 
and comparison with experiment. 

The SDI-deg results for Rsb are remarkably close to unity, ranging from 0.964 to 1.009. 
With FPD6, the results range from 0.826 to 0.989 — somewhat farther from unity, but again 
noticeably close. Remember that, in a very simple model for a vibrational nucleus, the ratio 
of Rsb would be zero. Clearly the shell model shows greater resistence for the intrinsic 
quadrupole moment Qo(S) to become small, as we might be led to believe from collective 
arguments. 

In Tables ITvT - IVIIH we show the results of M(Q)j and M(B)j^j_ 2 for all the nuclei with 
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FIG. 3: The same as Fig. [2] for 46 Ti. 

both interactions. We first look at M(Q)j for the surface delta interaction. For 44 Ti, 46 Ti, 
and 48 Ti, the values are slightly larger but remarkably close to unity, even for very high 
spins, e.g., I = 12, 16, and 18 for A = 44,46, and 48, respectively. In 48 Cr the results are 
not so close beyond 1 = 6, the values being 0.859, 0.766, and 0.779 for 1 = 8, 10, and 12±. 
Strangely, for 50 Cr the results are better up to / = 10. All in all, though, we are very far 
away from the simple vibrational limit of zero. 

If we look at M(-B)/_,.j_2 with SDI-deg, the results up to the transition 8^6 are all 
greater than 0.8 and less than unity for all nuclei considered. Beyond that, there are some 
lower results that may be due to band crossing, e.g., the value of M(£>) 12l _^ 10l in 46 Ti is 
0.0044; however, M(B) 12a ^i 0l is 0.6654. 
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FIG. 4: The same as Fig. [2] for 48 Ti. 

Results for the more realistic FPD6 interaction including single-particle energies are also 
shown in Tables lIVHVIIIl The results are over all not as close to unity as with SDI-deg. 
Still, one gets some substantial static quadrupole moments. Sometimes, the lowest state of 
a given angular momentum does not belong to a K = band, e.g., the I = 6f state in 48 Ti 
or the 10^ state in 50 Cr. The sign of the static quadrupole moment is opposite to what one 
would get assuming K = 0. In 48 Ti, the second 6 + state fits better into this profile. The 



near degeneracies of the 6^ and 69 states in Ti has been discussed previously [21]. 



In previous works on 50 Cr 22j, |23|, it was noted that the first 10 + state did not belong to 
the K = ground state band. Indeed, if one must choose a K value, it would seem K = 10 is 
the best for the 10^ state. This is supported by the fact that the static quadrupole moment 
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FIG. 5: The same as Fig. EI for 48 Cr. 

is large and positive, while the K = static moments are negative. Moreover, if it were 
strictly K = 10 and the 8^ were strictly K — 0, the B(E2) would be strongly inhibited. 
The small value of M(_B) 10l ^s — 0.2678 somewhat supports this. The decay lOj" — » 8^ is 
stronger, with Af(S)io 3 _» 8 = 0.6599. 

Let us briefly discuss Figs. [2H6] corresponding to 44 Ti, 46 Ti, 48 Ti, 48 Cr, and 50 Cr, respec- 
tively. Although not perfect, the FPD6 interaction in a full fp space yields a pretty good 
agreement for the energy levels of all nuclei here considered. The SDI interaction, for which 
the 0-2 splitting is fitted to experiment, gives a more spread out spectrum. It is closer to an 
/(/ + 1) spectrum than results with FPD6, but still significantly different. The spreading 
of the spectrum with SDI-deg is mainly due to the fact that there are no single-particle 
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FIG. 6: The same as Fig. [2] for 50 Cr. 

splittings in this model, i.e., e/ 7/2 = e/ 5/2 = e P3/2 = e Pl/2 . 

Just to give some numbers, in 48 Ti the experimental energy of the 12^ state is 8.09 MeV, 
FPD6 gives 8.31 MeV, and SDI-deg gives 15.98 MeV. Using the simple rotational model and 
fitting the 0-2 splitting to experiment, the 12^ state would, with an /(/ + 1) spectrum, be 
at 25.57 MeV. 

IV. THE HARMONIC VIBRATOR 

In the harmonic vibrator model, the nuclear shape oscillates between oblate and prolate. 
One gets equally spaced spectra, i.e., the ground state has angular momentum / = 0, the 
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TABLE IV: Ratios of intrinsic quadrupole moments M(Q)j and M(-B)j_»/_2 (see text) for 44 Ti 
using both SDI-deg and FPD6 interactions. 







SDI-deg 




FPD6 


I -► 


1-2 


M(Q)j 


M(S)i_j_ 2 


M(Q)j 


Af(S)j_j_ 2 


2 


-► 


1 


1 


1 


1 


4 


-> 2 


1.0092 


0.9741 


1.0534 


0.9764 


C 


-> 4 


1.0176 


0.9096 


1.1047 


0.8554 


8 - 


-> 6 


1.0276 


0.8241 


0.8498 


0.6971 


10 


-» 8 


1.0406 


0.7239 


0.7800 


0.6802 


12 - 


-► 10 


1.0556 


0.5498 


0.8472 


0.5130 



TABLE V: The same as Table ED for 46 Ti. 



SDI-deg FPD6 



/ -» 


1-2 


M(Q)j 




M(Q)j 


M(B):^< 


2 - 


-> 


1 


1 


1 


1 


4 


-> 2 


1.0159 


0.9878 


1.0349 


0.9975 


6 


-> 4 


1.0328 


0.9351 


0.9483 


0.9398 


8 - 


-► 6 


1.0295 


0.8152 


0.9522 


0.8626 


10 


8 


1.0270 


0.7562 


0.9362 


0.7436 


12i 


-> 10 


1.0131 


0.0044 


0.6763 


0.4493 


122 


-> 10 


1.0658 


0.6654 


0.2357 


0.2245 


14- 


-+ 12 x 


1.0331 


0.6803 


0.6724 


0.3804 


14- 


-+ 122 


// 


0.0007 


// 


0.2709 


16 - 


-> 14 


1.0480 


0.5212 


0.8911 


0.0644 



first excited state has angular momentum 1 = 2 and energy E(2). At 2E(2) there are states 
with I = 0, 2, and 4; at 3E(2) there are states with I = 0,2, 3, 4, and 6; etc. 

The selection rules and B(E2) , s relations are given in Bohr and Mottelson vol. 2, page 
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TABLE VI: The same as Table ED for 48 Ti. 







SDI-deg 




FPD6 


I -» 


1-2 


M(0)t 




M(Q)j 


M(B)j_j_ 2 


2 - 


-> 


1 


1 


1 


1 


4 


-> 2 


1.0172 


1.0015 


0.8654 


1.0310 


Oi 


-> 4 


1.0276 


0.9884 


-0.8738 


0.5016 


62 


-> 4 


0.2782 


0309 


0.8041 


0.7124 


8 - 


-> 61 


1.0319 


0.9598 


0.6134 


0.5768 




- 




V 


0.0016 


// 


0.3832 


10 


-» 8 


1.0340 


0.9193 


0.8153 


0.5682 


12 - 


-> 10 


1.0347 


0.8589 


0.5462 


0.3616 


14 - 


-> 12 


1.0392 


0.7821 


0.8447 


0.1491 


16 - 


14 


1.0519 


0.6657 


0.9128 


0.4386 


18 - 


-> 16 


1.0614 


0.5062 


1.0951 


0.0646 



349 24j. The equations for a transition are 



Cn-l^n-l 



B(EX; n x CJn -> nA-iCn-i4-i) = n x B(EX, n x = 1 -> n A = 0) , (9) 



i.e., there is stimulated emission — the more quanta there are, the bigger the B(EX). In 
the above, n> is the number of vibrational quanta, I n is the angular momentum and ( n 
stands for any additional quantum numbers. Bohr and Mottelson also give the selection 
rule An a = ±1, which implies that static quadrupole moments vanish, consistent with the 



opening sentence in this section 



24] 



Since we are considering transitions from a state with the maximum angular momentum 
(imax) of all degenerate states with n\ quanta, the final state for E2 transitions must have 
I = Imax — 2, so there is only one term in the left-hand side: J max — 2. Note that n\ = / max /2 

B(E2) M . 2 = I -B(E2) 2 ^ . (10) 



2 



Hence, we obtain 



1/2 
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TABLE VII: The same as Table ED for 48 Cr. 



SDI-deg FPD6 



I -> 


1-2 


M(Q)j 


M(B)j^j_ 2 


M(Q)j 


M(5)j_j_: 


2 - 


-> 


1 


1 


1 


1 


4 


-> 2 


0.9797 


0.9998 


1.0088 


0.9887 


6 


4 


0.9028 


0.9851 


0.9679 


0.9600 


8 - 


-► 6 


0.8586 


0.9592 


0.9370 


0.9101 


10 


-> 8 


0.7660 


0.9121 


0.7703 


0.8028 


12i 


-> 10 


0.7788 


0.8839 


0.1381 


0.5277 


12 2 


-► 10 


1.3175 


0.0985 


0.5786 


0.5102 


14! ■ 


-» 12i 


1.3722 


0.0790 


0.1675 


0.5024 


14 2 ■ 


-» 12 x 


0.7588 


0.8642 


0.6645 


0.1394 


16i ■ 


14i 


1.4145 


0.4311 


0.1537 


0.3536 


16 2 ■ 


-14i 


0.7379 


0.1255 


0.6231 


0.1401 


16 2 ■ 


14a 


// 


0.6032 


// 


0.5163 


18! ■ 


-+ 16i 


1.1348 


0.0000 


0.5938 


0.0577 


18! - 


-► 16 2 


// 


0.0000 


// 


0.4369 


18 2 ■ 


16i 


0.8025 


0.0000 


0.7394 


0.0157 


18 2 ■ 


-» 16 2 


// 


0.5765 


// 


0.1430 


20 ■ 


-» 18 


1.0753 


0.4473 


0.8097 


0.1183 



Note that M(5)/_ > /_2 increases steadily with J, which is not the case with SDI-deg 
or FPD6. Some values of M vih (B) I ^ I _ 2 for I = 0,2,4,6,8, and 10 are, respectively, 
1, 1.183, 1.381, 1.558, 1.719, and 1.867. For large values of /, we have 

M(B)^ 7 _ a - I -I J , (12) 

still growing steadily with I. 

As mentioned before, the vibrational prediction for static quadrupole moments is 
M v fo(Q)i = 0, which is certainly not the case with SDI-deg or FPD6. 
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TABLE VIII: The same as Table ED for 50 Cr. 



SDI-deg FPD6 



I -> 


1-2 


M(Q)j 


M(B)j^j_ 2 


M(Q)j 


M(5)j_j_: 


2 - 


-> 


1 


1 


1 


1 


4 - 


-> 2 


0.9881 


1.0027 


1.0014 


1.0074 


6 - 


-> 4 


0.9631 


0.9954 


0.5587 


0.8785 


8 - 


6 


0.9332 


0.9807 


0.6464 


0.8154 


10i 


-> 8 


0.9064 


0.9557 


-0.9601 


0.2678 


10 2 


-► 8 


1.0923 


0.0094 


0.2300 


0.6599 


12i - 


-+ 10i 


0.8857 


0.9751 


-0.3819 


0.3115 


12i - 


-+ io 2 




0.0959 




0.2502 


12 2 - 


-+ 10i 


1.1533 


0.0058 


0.5044 


0.1502 


12 2 - 


-+ io 2 


// 


0.8039 




0.5064 


14i - 


12i 


0.8482 


0.8672 


-0.2057 


0.3662 


14i - 


12 2 




0.1310 




0.1013 


14 2 - 


-+ 12i 


1.1224 


0.0131 


0.3175 


0.0947 


14 2 - 


- 12 2 


// 


0.1404 


// 


0.4982 


16 - 


- 14i 


0.8832 


0.8584 


0.2218 


0.0943 


16 - 


14 2 


// 


0.0161 


// 


0.5892 


18 - 


16 


1.1256 


0.0266 


0.1888 


0.4213 


20 - 


-> 18 


1.1306 


0.6272 


0.5944 


0.0619 


22 - 


20 


1.0694 


0.4435 


0.9263 


0.1189 



V. CLOSING REMARKS 

This work is an extension of previous work by Robinson et al [lj], where it was noted 
that, for a large variety of nuclei, the simple rotational formula, if fitted to the experimental 
B(E2) from the ground state to the 2± state, could give a good result for the static quadruple 
moment of the 2^ state, i.e., in the notation of this work, Rsb is close to one. The one 
exception is 40 Ar. 
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In this work, we extend the calculations to higher energy and higher angular momentum. 
We use both a phenomenological interaction FPD6 with realistic single-particle splittings, 
and, as a counterpoint to Q ■ Q, we use the schematic surface delta interaction. We claim the 
latter shows also collective properties. As seen in Fig. [T], its spectrum is not rotational but 
it does have features that many nuclei possess — a collective appearance; and it should be 
noted that the nuclei we consider do not have rotational spectra either. We do not include 
single-particle splitings in SDI, and this is the main reason, rather than the SDI per se, that 
the spectrum, though less spread out than in the rotational model, is more spread out than 
with a realistic interaction or, indeed, experiment. The effects of single-particle splittings 
are taken care of in the realistic case. 

The results with SDI for M(Q)j (ratios of intrinsic quadrupole moments) is close to one 
for many nuclei and many angular momenta /. One only runs into trouble when one has 
near degeneracies like 6f and in 48 Ti, as well as lOt and lOt in 50 Cr. These degeneracies 



have been addressed previously by Zamick et al. 2l|, |22|, |23j]. And when one gets band 
crossings, the situation can get confused. 

But still, all in all, with SDI we get some remarkable agreements with the rotational 
formulas for ratios between static quadrupole moments and B(E2) , s, not only for angular 
momentum 1 = 2, but for higher / as well. The realistic interactions also yield similar 
results, although perhaps not quite as definitive as does the schematic SDI. 
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